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Abstract 

In this paper we make an attempt to give a consistent background and definitions suitable 
for the theory of integrable difference equations. We adapt a concept of recursion operator to 
difference equations and show that it generates an infinite sequence of symmetries and canonical 
conservation laws for a difference equation. Similar to the case of partial differential equations 
these canonical densities can serve as integrability conditions for difference equations. We have 
found the recursion operators for the Viallet and all ABS equations. 

1 Introduction 

The nature of integrability of partial differential equations is rather well understood [HE]- It can 
be related to the existence of Lax representations, Miura and Backlund transformations, bi-linear 
Hirota representations, multi-soliton solutions, the Painleve property, bi-Hamiltonian structures, 
recursion operators and infinite hierarchies of symmetries and conservation laws. Although there is 
not any commonly accepted definition for integrability, the patterns of equations possessing some of 
the above listed properties coincide. In the symmetry approach the existence of an infinite hierarchy 
of symmetries and conservation laws is taken as a definition of integrability and the most complete 
results on classification of integrable equations are based on this property O HI [5j [6j [H [9] . The 
symmetry approach proved to be effective in the problem of classification of integrable differential- 
difference equations [101 E] ■ 

The theory of integrable difference (or discrete) equations is a relatively recent, but very ac- 
tive area of research. It inherited some methods and tools from the theory of integrable partial 
differential equations, see for instance [1'2\ 113], The property of multidimensional consistency, or 
consistency around the cube [14] . has been proposed as the integrability criterion for quadrilateral 
difference equations [15, [TB]. All integrable affine-linear quadrilateral equations (satisfying a few 
symmetry conditions with respect to reflections of the three dimensional lattice) were classified by 
Adler, Bobenko and Suris [IT], see also [18]. They produced a list of seven equations (up to point 
fractional-linear transformations), which will be referred as the ABS equations. Unfortunately, the 
method based on the property of multi-dimensional consistency cannot be directly applied to a 
given equation on a lattice. Indeed, to check the consistency around the cube one needs to know 
the equations corresponding to the sides of the cube. 

Another characteristic property of an integrable equation is the vanishing of its algebraic entropy 
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[19J. It has been used by Viallet to identify the equation 

Q := aiUofiUifiUo^u^i + a 2 (u ,o'"i,o"o,i + "i,o"o,i"i,i + Wo,i«i,imo,q + "i,i"o,o"i,o) 

+O3(lto,0«l,0 + «0,l«l,l) + O4(«l,0W0,l + «0,0Wl,l) + 05^0,0^0,1 + "l,0«l,l) 

+06(^0,0 + «i,o + "0,1 + "i,l) + a 7 = 0, (1) 

where dj are free complex parameters [20] . By a point fractional- linear transformation equation 
flU with a generic choice of parameters can be reduced to Adler's equation, also referred as the Q4 
equation in the ABS classification. Therefore it is not surprising that the Viallet equation ([TJ has a 
hierarchy of symmetries [21] , All of the ABS equations can be obtained from the Viallet equation 
by a simple specialisation of parameters, it is a convenient way to represent the ABS equations. 
The obvious advantage of the Viallet approach is its applicability to a given equation. 

The aim of this paper is to develop a theory of difference equations, based on the study of 
symmetries, conservation laws and (formal) recursion operators. By integrability of a difference 
equation we understand the existence of an infinite hierarchy of its symmetries. In this paper 
we consider quadrilateral equations Q = 0, where Q is assumed to be an irreducible affine- linear 
polynomial over C. Symmetries, conservation laws, and actually all rational functions which make 
sense on solutions of the difference equation are elements of a difference field of fractions Tq defined 
in Section 2. The affine-linearity of the equation enables us to introduce the elimination map £ , 
which reduces any element of Tq to a rational function of dynamical variables. This map is quite 
a useful tool for computations and proofs, since two functions / and g of the lattice variables 
are equivalent modulo the difference equation if and only £(f) = £(g)- In Section 2 we give a 
self-contained set of basic definitions required in this paper. 

In Section 3 we give a definition of a recursion operator for difference equations. The recursion 
operator acting on a seed symmetry of the difference equation generates an infinite hierarchy of 
symmetries of that difference equation. We show that the residues of the powers of the recursion 
operator are (canonical) densities of conservation laws for the corresponding difference equation. 
We can generate an infinite sequence of canonical conservation laws (both the densities and fluxes). 
Similar to the differential case (see [3j HI [5] ) , these canonical conservation laws can be used as inte- 
grability conditions for a difference equation. They are the necessary conditions for the existence of 
a (formal) recursion operator and, ultimately, for an infinite hierarchy of symmetries of a difference 
equation. 

We present the recursion operators for the Viallet equation and all of the ABS equations in 
Section 4. These recursion operators generate infinite hierarchies of symmetries for every ABS 
equation. The corresponding canonical conservation laws provide us with an infinite sequence of 
densities and corresponding fluxes of increasing order. The very first term in this sequence was 
found by Hydon and Rasin for all of the ABS equations by a direct computation |22j. It has been 
shown recently that the Gardner method, originally developed for integrable partial differential 
equations, can be adapted to difference equations and it was applied for the production of infinitely 
many conservation laws for HI and some other equations, including non-symmetrical ones [23} 124], 

Our work was partially inspired by the results obtained by Levi and Yamilov. In [25J they 
started to develop a theory of integrable difference equations based on the study of their symme- 
tries. The goal of Levi and Yamilov theory is a classification of difference equations possessing 
symmetries. They made a good progress in a formalisation of the theory and they also discussed 
the issue of dynamical variables. Levi and Yamilov have found one integrability condition for dif- 
ference equations and it does coincide with the first element of our infinite sequence of canonical 
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conservation laws. Their approach is based solely on symmetries and does not use the concept of 
(formal) recursion operator. 

We would like to stress that our assumption, that the difference equation is an affine-linear 
polynomial and its symmetries are some rational functions, has been made mainly for the clarity 
of presentation, and it is suitable for the ABS equations. This assumption can be removed under 
some inessential technical conditions, which guarantee that the difference equation has a unique 
solution for an open set of initial data. Obviously all results in this paper are valid for local analytic 
functions of the lattice variables. 



2 Basic definitions 



In this section, we introduce the elimination map £ and define dynamical variables for quadrilat- 
eral affine-linear difference equations. To make the paper self-contained, we also give the basic 
definitions of symmetries and conservation laws of difference equations. 



2.1 Difference equations and dynamical variables 

Difference equations on Z 2 can be seen as a discrete analog of partial differential equations. Let 
us denote by u = u(n, m) a complex-valued function u : Z 2 h- > C where n and in are "independent 
variables" and u will play the role of a "dependent" variable in a difference equation. Instead of 
partial derivatives we have two commuting shift maps S and T defined as 

S : u i-)- uifi = u(n + 1, m), T : u h-> Uo,i = u(n, m + 1), S p T q : u h-> u ps = u(n + p,m + q). 

For uniformity of the notation, it is convenient to denote the "unshifted" function u as Uo,o- In the 
theory of difference equations we shall treat u P) q as variables. The set of all shifts of the variable u 
will be denoted by U = {u PtQ \ (p, q) £ 7?}. For a function / = f(u Pltqi , . . . , u PkSk ) of variables u P)q 
the action of the operators <S, T is defined as 

>5 T J '(f) = fi,j = f( u pi+i,qi+ji ■ ■ ■ } u Pk+i,qk+j) ' 

A quadrilateral difference equation can be defined as 

Q(uo,o, "l,o, "0,1,1*1,1) = , (2) 

where Q(uo,o> "1,0) "0,1, "1,1) is an irreducible polynomial of the "dependent variable" u = uo,o and 
its shifts. Irreducibility means that Q cannot be factorised and presented as a product of two 
polynomials. It is assumed that equation J2|) is valid at every point (n, m) £ 1? and thus (|2|) 
represents the infinite set of equations 

Qp,q = Q(up, q ,Up +liq ,Up tq+ i,u p+liq+1 ) = 0, (p, q) G I? . (3) 

In this paper we shall consider only polynomials Q. This restriction is not very essential for our 
construction, but allows us to make the presentation more transparent. Moreover, for quadrilateral 
equations we shall assume that Q is an irreducible affine-linear polynomial which depends non- 
trivially on all variables, i.e. 

QQ QlQ 

7^0, = 0, i,j £ {0, I}, Q G C[«o,0,"l,0,"0,l,«l,l]- 
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This is true for all of the ABS equations (except the trivial case of the HI equation when a = /?). 

Let C[U] be the ring of polynomials of the variables U, where it is assumed that every polynomial 
has a finite number of terms. Maps S and T are automorphisms of C[U] and thus C[U] is a difference 
ring. We denote Jq = ({Q p , q }) the ideal generated by the difference equation and all its shifts ©. 
Since Q is an irreducible affine linear polynomial, the ideal Jq is prime and radical. For any a G Jq 
we have S(a) G Jq and T(a) G Jq. Thus Jq is a difference ideal and C[U]/Jq is a difference 
quotient ring. 

The fields of rational functions of variables 

U a = {u nfi | n G Z}, U t = {u ,„ | n G Z}, f/ = ^ s U I7 t . 

are denoted as 

F s = C(U S ), F t = C(U t ), F = C(U ) 

respectively. Rational functions of variables u Pt g that can be restricted to the solutions of the 
difference equation (i.e. the denominators do not vanish on the solutions) form a field 

F Q = {[a]/[b]\a,beC[U], b$J Q }, 

where [a] denotes the class of equivalent polynomials (two polynomials f,g€ C[U] are equivalent, 
denoted by / = g, if / — g G Jq). For a,b,c,d G C[U], b,d Jq, the two rational functions 
a/b and c/d represent the same element of J-q if ad — be G Jq. Such rational functions we shall 
call equivalent and denote a/b = c/d. The values of equivalent rational functions coincide for all 
solutions of the difference equation. 

In general, it is quite a tedious computational problem to work with polynomials and rational 
functions modulo an ideal. In our case the ideal Jq is generated by affine-linear polynomials Q p , q 
and we can uniquely resolve equation Q = with respect to each variable 

^0,0 = F{uifl, « 0) 1 ), 1*1,0 = G^mq.o, fo.i, /.n, 

"0,1 = #0*0,0) ^1,0, «l,l)> "1,1 = MOiQ ) o,Ml,0,«0,l)- 

Since Q is an affine linear polynomial, functions F, G, H and M are rational functions of their 
arguments. Equations dH) enable us recursively and uniquely to express any variable u Ptq in terms 
of the variables Uq = U s U Ut- 

Definition 1. For elements of U the elimination map £ : U i— >■ C(Uq) is defined recursively: 

Vp G Z, £0*o,p) = «o, P , ^( u p,o) = «p,o, 

if p > 0,g > 0, £(u P) g) = M(£ , (n p _i >g _i),£'(np i q_i),^(up_i i(? )) , 

ifp<0,g>0, £(u P ,q) = H(£(u Pl g-i),£ (u p+1)q - 1 ),£(u p+1<q )) , (5) 

if p > 0, q < 0, £(u P ,g) = G(£(up-i )q ),£(up-i >q+ x),£(itp )q+ i)) , 

if p < 0, g < 0, £(up !q ) = F(£{u p+ i >q ),£{u PA+ i),£(u p+ ^ q+1 )) . 

For polynomials f(u Pltqi , . . . ,u Pktqk ) G C[U] the elimination map £ : C[U] i-> C(J7o) is defined as 

)) G C(C/ ). 

For rational functions a/b, a, b G C[£7], b ^ Jq the elimination map £ is defined as 

£ : a/b^ £(a)/£(b). 
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Variables Uq we shall call the dynamical variables, it is easy to see that £(u Ptq ) is a rational 
function of \p\ + |g| + 1 dynamical variables 

S(u Pt g) £ C({u nfi ,u 0jm |0 < \n-p\ < \p\, < \m - q\ < \q\}). (6) 

The choice of the set of dynamical variables is not unique. Another popular choice is the so-called 
"staircase" set {u ntn ,u n ^ n+ i \ n € Z} which will not be used in this paper. 

The elimination map £ : C[U] \— > C(Uq) is a difference ring homomorphism, its kernel is Jq, and 
its image is a difference ring isomorphic to the factor ring C[U]/Jq. The field C(Uq) is a difference 
field with automorphisms £ o S and £ o T ■ The elimination map £ : Tq h- > C(Uq) is a difference 
field isomorphism. It is a useful tool to establish whether two rational functions /, g of variables U 
are equivalent (i.e. represent the same element of Tq): f = g if and only if £(f) = £(g). Later on 
we shall use logarithms of rational functions. The action of £ on ln(a), a £ J-q, a Jq is defined 
as £(ln(a)) = ln(£(a)) and thus we shall say that ln(a) = ln(6) if and only if ln(£(a)) = ln(£(b)). 

To illustrate the above construction, let us consider the difference equation Q = with 

Q = ("o,o -«i,i)(«i,o - u ,i) - at + P, a,/3eC, (7) 

which is known as the HI equation in the ABS classification or the discrete potential KdV equation. 
First we should assume that a ^ (3, otherwise the polynomial Q is reducible. It follows from 
Definition [1] that 

a — f3 

£{u^i) = M(u ,o,ui,o, «o,i) = ^o,o 



«1,0 - "0,1 

s Uq l"0,2 - "0,0"o 1 ~ ( a ~ P) U 1,0 + "0,0"0,l"l,0 ~ "0,l"0,2"l,0 
£{Ul,2) = ] a ' ■ 

P - a - UQflU 0; i + -Uo,1^0,2 + "o,o"i,o - "0,2"l,0 

Two rational functions 

, _ "o,o - "1,1 g_ a- P 

"i,o ~ u-i,o ' ("i,o ~ "-i,o)(«i,o - "o,i) 

are equivalent / = g since 

uofl ~ "1,1 a- P Q 



"i,o-"-i,o ("i,o -"-i,o)(«i,o -«o,i) («i,o -«-i,o)("i,o -"o,i) 
and it is easy to check that £(/) = £(g) = g- 

2.2 Symmetries and conservation laws of difference equations 

Definition 2. Let Q = be a difference equation. Then K € Tq is called a generator of an 
infinitesimal symmetry ( or simply, a symmetry ) of the difference equation if 

D Q (K)^0. (8) 

Here Dq is the Frechet derivative of Q defined as 

D Q = Y J Qu i , j s i T j , Q^ = ^-. (9) 



5 



The sum in ([9]) is finite for a given Q. For a quadrilateral equation it has only four terms. Here 
Dq(K) is equal to zero as element of J-q, i.e. Dq{K) C Jq. The way to check it is to apply the 
elimination map: thus, K is a symmetry if £(Dq(K)) = 0. 

If K is a symmetry and u = u(n, m) is a solution of a difference equation ([2]), then there is an 
infinitesimal transformation of solution u: 



u = u + 



satisfying 

Q(uo,o,Ui t0 ,u ,i,u lt i) = 0(e 2 ). 

If the difference equation Q = admits symmetries, then they form a Lie algebra. Indeed, with 
a symmetry K G J-q we can associate an evolutionary derivation of J-q (or a vector field on ^q): 

Xk= £ ^asb' 

where = S p T q (K). A derivation (a vector field) Xr- is called evolutionary if it commutes with 
the shift automorphisms: SXk = XkS, TXk = XkT. It follows from (jHJ) that for any a £ Jq we 
have Jjf(o) G Jq and thus the evolutionary derivation Xk is defined correctly on the difference 
field of fractions J-q. 

If the difference equation admits two symmetries F and G, then the commutator of the corre- 
sponding evolutionary derivations is an evolutionary derivation of J-q 

XfXq — XqXf = Xh, 

where H = [F, G], with [F, G] denoting the Lie bracket 

[F, G] = X F (G) - X G {F) = D G (F) - D F {G) G Tq. (10) 

It is easy to see that H is a symmetry, indeed 

D Q (H) = X H (Q) = X F X G {Q) - X G X F {Q) = X F D Q (G) - X g Dq(F) = 0. 

Symmetries form a linear space over C. The bracket (jlOp is linear with respect to each argument, it 
is skew-symmetric and satisfies the Jacobi identity. The Lie algebra of symmetries of the difference 
equation Q = will be denoted as 21q. Existence of an infinite dimensional Lie algebra 21q is a 
characteristic property of integrable equations and can be taken as a definition of integrability. 

We would like to emphasise that a symmetry K is an element of J-q, so it can be represented 
by a function 8{K) G J-q of variables Uq. As an element of Fq, a symmetry K depends on a finite 
set of variables {ii n ,o I -^l < n < N2} and {uo,m I Mi < n < M2}, and thus can be characterised by 
four integers (N\, N2, Mi, M2). It is straightforward to check that, if 

K = K(un u o, ■ ■ ■ ,un 2 ,o,uo,m 1 , ■ ■ ■ ,Uo,M 2 ) 

is a symmetry of a quadrilateral equation, then it is a sum of two functions 

K = K s (u Nlfi , ... , u N2i o) + Ktiuo^Mx, ■ ■ ■ , uo,M 2 ) ■ 
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In particular, for the ABS equations it is known |21|, [26] that, any "five point symmetry" K = 
K(u ,o, U-1,0, u _i, «i )0 , u ,i) is a sum of two symmetries K = K s (u-i )0 , u ,o, «l,o)+-^t(«o «o,o, «o,l)- 
Thus, the problem of the symmetry description is likely to be reduced to the study of the symme- 
tries depending on "one-dimensional" set of variables U s or Ut- Moreover, all ABS equations and 
the Viallet equation possess infinite hierarchy of such symmetries. 

Let / = /(itjVi,o> • • • !^at 2 ,o) ^ T s and df/duN lt o ^ 0> df /duN 2) o ^ 0, then the order of / is 
defined as ord s (/) = (Ni,Nz). For example, 

= (ui,o - u-1,0)" 1 (11) 

is the generator of a symmetry for equation (|7|) and oid s (K^) = (—1,1). Similarly, the order of 
9 = g(uo,Mi > • • • > n o,M 2 ) € J" t is defined as ord t (#) = (M 1 ,M 2 ). 

Definition 3. 1. A pair of functions p,a G J-q is called a conservation law for the difference 
equation (0), if 

(T-l)(p) = (S-l)(a). (12) 

Functions p and a will be referred to as the density and the flux of the conservation law and 
1 denotes the identity map. 

2. A conservation law is called trivial, if functions p and a are components of a (difference) 
gradient of some element H G Tq, i.e. 

p=(S-l)(H), a = (T-l)(H). 

3. A non-constant element C € J-q is called s-constant, (respectively t-constant) if S{C) = C 
(respectively T(C) = C). Elements of J-q which are invariant with respect to the both shifts 
are called constants of Tq. 

Existence of non-trivial conservation laws and/or s or t-constants for a difference equation is a 
rather exceptional property. The ABS equations (and many other integrable difference equations) 
possess infinitely many nontrivial conservation laws. Similar to the differential case [27J, the ex- 
istence of s or t-constants suggests that the difference equation is linearizable. For instance the 
linear difference equation with Q = uo,o — ^1,0 — u o,i + u \,i has a t-constant f(ui t o — Uo,o) an d a 
s-constant 3(1*0,1 — ^o,o)> where /, g are arbitrary functions. 

Definition 4. Let Im(5 — 1) and Im(T— 1) denote the images of the maps S — 1 : J-q h4 J-q and 

T — 1 : Tq 1— >■ Tq respectively. 

1. Two elements p, g £ Tq are s-equivalent (p = s g) if p— g € Im(5 — 1). Elements of the factor 
space Tq/\vh(S — 1) we shall call s-functionals, or densities. 

2. Two elements a, c € Tq are t-equivalent (a =t if o~ — s G Im(T — 1). Elements of the 
factor space Tq/Ivh^T — 1) we shall call t-functionals, or fluxes. 

In general the problem to establish whether two elements of Tq are s- or t-equivalent is diffi- 
cult. This problem can be easily solved for densities belonging to T s (or fluxes belonging to Tt). 
Moreover, in the next Section we will show that an integrable quadrilateral equation (such as the 
the Viallet equation) possesses an infinite hierarchies of canonical conservation laws such that their 
densities are elements of T s . 

There is a criterion to determine whether two elements of T s are equivalent or not. This is 
based on the notion of variational derivative (Euler's operator). 
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Definition 5. Let / G 7 S has order (A^i,A^), then the variational derivative 5 S of f is defined as 



N 2 /at 

df 



k=Ni 



Two elements p, g £ T s are equivalent p = s g if and only if 5 s (p) = S s (g). In particularly, if 
density p € T s is trivial, then 8 s (p) = 0. The order of a density p G T s is defined as ord<5 s (p) = 
N2 — N±, where (N\,N2) = ord(5 s (p)). Equivalent densities have the same order. For example the 
densities of conservation laws for the HI equation (|7|) 

po = ln(ui,o - U-1,0) and pi 



(ui,o - u_i j0 )(w2,o - Uo,o) ' 
are of orders 4 and 6 respectively. 



3 Recursion operator and canonical conservation laws 

In this section we give definitions of (s- and t- ) difference and pseudo-difference operators. We 
adapt the concept of recursion operators to difference equations and give the criteria for a pseudo- 
difference operator to be a recursion operator of a quadrilateral difference equation. Using Laurent 
formal series and their residues we derive infinitely many integrability conditions for quadrilateral 
difference equations. 



3.1 Recursion operators of difference equations 

In order to proceed we need to define s- and t-pseudo-difference operators. 

Definition 6. (1) A s-difference operator B of order N with coefficients in Tq is a finite sum of 
the form 

B = b N S N + ^icS^- 1 + • • • + b M S M , b k € Tq, M < N, N, M e Z. (14) 

(2) An operator L is called a s-pseudo- difference operator if it is either a finite sum of terms 

L = Y,B n oC~ l (15) 

(where B n and C n are s-difference operators with coefficients in Tq and o denotes their composition) 
or a finite composition of s-pseudo- difference operators of the form (T5\). 

In the above definition we used the shift map S, and thus we have defined the s-difference, 
s-pseudo-difference operators. Using the shift map T instead, one can define the t-difference and 
pseudo-difference operators. We shall omit s- (or t-) and simply say difference or pseudo-difference 
operators if it is clear from the context. 

There is a natural action of difference operators on all elements of Tq but the action of pseudo- 
difference operators is not generally defined. For a pseudo-difference operator L (|15p it can be 
defined on those elements of a» G Tq that belong to the intersection of the image spaces of the 
difference operators C n , so that the action of operators C~ 1 (ai) is defined, i.e. C~ 1 (ai) € Tq. A 
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linear space that is a C-span of all such elements a% is a domain of the pseudo-difference operator 
denoted by Dom(L). 

By a recursion operator 9\ of a difference equation Q = we shall understand a s-pseudo- 
difference operator such that 

m : Dom(9l) n 21q H- 21q, 

where 21q is the linear space of symmetries of this difference equation. In other words, if the 
action of on a symmetry K 6 Tq is defined, i.e. dK{K) G Tq, then 9^(ET) is a symmetry of 
the same difference equation. Similarly, we define a t-pseudo-difference recursion operator *H : 
Dom(^H) PI 21q h-> 21q. The operator of multiplication by a constant is a trivial recursion operator. 
In what follows we shall assume that 9\, 9\ are nontrivial. 

Theorem 1. Let Q(uQ t Q, u^o, uo,l> u i,i) = 6e a difference equation. 

(i) if i/iere exist two s-pseudo- differential operators 9^ and such that 

Dqo^ = ^oDq, (16) 

t/ien 9\ is a recursion operator of the difference equation. 

(ii) Relation Mb}) is valid if and only if 

T(9V)-<R= [Qoft,®- 1 } , (17) 

where & = (Q Ul 1 <S + Q UQ i) 1 ° (Qui <S + Qu 00 ), and the operator *P satisfies 

= (Q Ulfi S + Q U0>0 ) o m o (Q U10 S + Q UQ , o y l - (18) 

Proof, (i) If K is a symmetry and SK(ir) G Fq it follows from CED that Dq2K(K) = <#D Q (K) = 
and thus 9\(K) is a symmetry. 

To prove (ii), we represent Dq in the factorised form 

Dq = Q Ultl S T + Q UoA T + Q uifi s + Qu ,o = (Qu lt iS + <5« ,i) ( 7 " + $ ) ' 

where 

$ = (Q« M S + Q^i) -1 ° (Q« ll0 5 + Q«o, ). (19) 

Then (|16|) can be rewritten as 

T(<K)oT + $o<K= (Qn 1 , 1 5 + Q Moa )- 1 oq}o(Q Ul l 5 + Q U(U )o(r+$) , (20) 
which is equivalent to 

T(3t) = (Q« 1 ,x5 + g« ,ir 1 oqJo(Q«i ) i5 + Q w , 1 ) (21) 

and 

$ o <H = (Q U1>1 5 + Q^J" 1 o <p o (Q U11< S + Q U0 l ) o (22) 
Using CED we can rewrite J22J in the form (USD . It follows from (HQ and ([22]) that 

$ o JK = r(9t) o $ (23) 



The latter can be rewritten in the form (TTj 

In the same way one can prove the similar statement in the t-direction. 



9 



Theorem 2. Let Q(^o,o ; u i,0i u o,i, u i,i) = be a difference equation. 

(i) If there exist two t-pseudo- differential operators and *p such that 

D Q o<h = yoD Q , (24) 

then 9\ is a recursion operator of the difference equation. 

(ii) Relation (24\) is valid if and only if 

- & = [* o ft, f - 1 ] (25) 
where ^ = (Q Ul X T + Q Ul ) _1 ° (Qu i"T + Qu )> and the operator can be written as 

*P = (Qu , 1 T + Qu 0fi ) {Qu 0:1 T+ Quofi)" 1 - 

As in the case of partial differential equations, we show that if 9\ is a recursion operator, so is 
9\ n for all n G Z in the following statement. 

Corollary 1. 1. Under the conditions of (i) in Theorem [7J the pseudo- difference operator 1H 
satisfies the following equations 

T(JK n ) - <K n = [$ o , S- 1 ], n G Z, (26) 

2. Under the conditions of (i) in T/ieorem t/ie pseudo- difference operator D\ satisfies equations 

S(9T) -<K n = o^,*" 1 ] nGZ, (27) 

Proof. For the first part, it follows from (fl~6|) that Dq o $H n = o Dq,u G Z. Thus, we can 
apply Theorem Q] to 5K n and *}} n to produce (|26p from (|17p . The proof of the second part of the 
Corollary is similar. ■ 

In all our definitions and statements there is an obvious symmetry between "T" and "<S" objects 
(such as symmetries, pseudo-differential operators, symplectic and Hamiltonian operators). In what 
follows we shall concentrate on "<S" objects, but formulate some statements for both if required. 

3.2 Canonical conservation laws and integrability conditions 

Definition 7. A formal (Lauren^) series of order N is defined as a formal semi-infinite sum 

A = a N S N + a N - 1 S N - 1 + ---+a 1 S + a + a_ 1 S- 1 + --- , a N ^ 0, a k G Tq, N G Z. (28) 

Laurent formal series form a skew-field. Sums and products (compositions) of formal series 
are formal series. The product is associative, but not commutative. For any formal series A 
there exists a formal series A -1 such that A o A' 1 = A" 1 o A = 1. In order to find the first n 
coefficients of A -1 one needs to know exactly the first n coefficients of A. The same holds for the 

One can also define a formal Taylor series of order — N as a semi-infinite sum 

C = c- N S~ N + ci_jv5 1_iv H + c— i<S — 1 + co + ci«S H , cat / 0, c fc e „Fq, NeZ 

to develop a similar theory based on the Taylor series (or both). But in the case when the recursion operator is a 
composition of a Hamiltonian and symplectic operators it can be shown that the both approaches give equivalent 
results. Thus in applications to the ABS and Viallet equations we restrict ourselves to Laurent formal series only. 
Study of the both (Laurent and Taylor) series make sense for asymmetric difference equations. 



10 



composition of two formal series. For instance, for two Laurent formal series A defined by (|28p . 
and B = b M S M + & M -i«S M-1 • • we have 

oo N+M-k 

AoB = J2 c N+ M-kS N+M ~\ c k = a N ^ q S N -o(b q+k ^ N ) G T Q . 

k=0 q=0 

Any pseudo-difference operator B can be uniquely represented by a Laurent formal series B^. 
For example for B = (aS + 6) _1 we have 

B L = a^S- 1 + a- 2 S~ 2 + a- 3 ^ 3 + • • • , 

where the coefficients a k G Tq can be found recursively: 

Definition 8. Let ^4^ denote the Laurent series representations of a pseudo-differential operator 
A. Then oidA^ is called the Laurent order of A. 

Now we introduce the residue of a formal series, as well as a discrete version of Adler's theorem 
|28| , which will be used consequently to produce integrability conditions for quadrilateral difference 
equations. 

Definition 9. Let A defined by (28\) be a formal series of order N . The residue ies(A) and the 
logarithmic residue resln(^4) are defined as 

res(A) = ao, resln(A) = ln(ajv) . 

We shall need the following difference analog of the Adler Theorem [28J : 

Theorem 3. Let A = a N S N + ajv-iS^- 1 • • • and B = b M S M + b M -iS M ~ l ■■■be two Laurent 
formal series of order N and M respectively. Then 

res[A,B] = (S-l)(a(A,B)), 

where 

N n M n 

a(A,B) =Y,Y,S- k (a~n)S n ~ k (bn) - £ J>- fc (&-„)<S n ~ fc K) G T Q . 

n=l k=l n=l k=l 

Proof. Only commutators of the form [a p S' p , b^ p S~ p ], p G Z contribute to the residue. As- 
suming p > 0, we get 

v 

[op5P,6_ p 5-P] = a p S p (b^ p )-S- p (a p )b^ p = (S p -l)(S- p (a p )b„ p ) = (5-1) £ S p - k (S- p (a p )b^ p )M 

k=i 

For first order formal series A the first n coefficients of the series completely define the first n 
residues po = res In A, p\ = resA, . . . , p n ~ l = resA n_1 . 

In the theory of pseudo-differential formal series (see for example [31 G3 E]) for a series A = 
aD% + bD^~ l + ■ ■ • of order N it is always possible to find the iV-th root, i.e. a formal series 
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B = An = aiv D x + (jjaN 1 b — ^j^-a^ 1 D x (a)j + • • • such that B N = A and, consequently, to 

parameterize the coefficients of A by a sequence of canonical residues res In B, resB, res-B 2 , .... In 
the case of formal difference series in general it is not possible to find the iV-th root of a series 
A = a^S N + aj\f^iS N ~ 1 + • • • , i.e. a first order series B = a±S + ao + ■ • • with the coefficients in 
J-q, such that B = A. Indeed, computations show that to determine the coefficient a k one needs 
to solve equations of the form 

N-i 

Y[ S i (a l ) = a N , a k + S(a k ) + ---+S N ~ 1 (a k )=g k , k = 0, -1, -2, . . . 

j=0 

for some g k G Tq. Obviously, a solution a k G Tq does not exist if g k ^ Im(l + S + • • • 5 7V_1 ). 

Theorem A. If a difference equation possesses a recursion operator £R, oic\l{0\) = N > 0, then it 
has infinitely many canonical conservation laws 

(T- l)pnN = (<S - l)cr nN , n = 0,1,2,... 

with canonical conserved densities 

po = reslnlKi, p n N = res£R£, n > 0. (29) 

Proof. Equation fj26[) is obviously satisfied if we replace and OK by the corresponding Laurent 
formal series &l and OKl- For n > 1, we compute the residue of equation (|26p to arrive at 

T(res^)-res^ = res([<£ L o^, ^}) = (S - l)a(^ L oOKl,^ 1 ). (30) 
It follows from Theorem [3] that 

a nN = a{^ L odVl,^l 1 ), a n eF Q . (31) 

For n = , we rewrite fj23[) as 

$ o m o o or 1 = t(jh) o ?r 1 . 

and compute its logarithmic residue. This leads to 

N-l 

T(resln<K L ) - resln<H L = (S - 1) ^ S^lnao), (32) 

fc=0 

where ao is the first coefficient in the expansion = ao + ai5 _1 + ■ ■ ■ , (see ([37|) ). ■ 
In the above Theorem, densities and fluxes are counted modulo N, where N is the Laurent 

order of the recursion operator. 

If a recursion operator is known, then Theorem [J] gives us a completely algorithmic way to find 

explicitly a sequence of conservation laws, including both the densities p k and the corresponding 

fluxes a k . The residues of powers a formal series are easy to compute. For instance, if 

OK = ri S + r + r-iS- 1 + r_ 2 S" 2 + r_ 3 S~ 3 + 

then 

reslnlH = lnri, resOK = ro, res9^ 2 = <S -1 (ri)r_i + r 2 . + rnS -1 (r_i), 
res<K 3 = 5- 2 (ri)5~ 1 (ri)r_ 2 + 5" 1 (r )5" 1 (ri)r_i + 25" 1 (ri)r_ir (33) 
+S- 1 (r 1 )r 1 S(r_ 2 ) + r 3 + 2r r 1 S(r_ 2 ) + r 1 «S(r_i)«S(r ) + ^(n^V-a). 
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Proposition 1. If a recursion operator £ft is represented by a first order formal series 9\l = 
r\S + ro + r_i5 _1 + • • • , then 

(i) (T - 1) (In n) = (5 - l)*" 1 f In ^ V (34) 



(ii) (T-l)(r ) = (5-l)5- 1 (r 1 F), (35) 

(iii) (T-l)(r_ 1 5- 1 (r 1 ) + r 2 + r 1 5(r_ 1 )) = (5-l)(a 2 ), (36) 

a 2 = 5- x (n F) {S-\r ) + r - S~ 2 (n F)} - (1 + S" 1 ) (n (n F)) , 

and F, G denote 

Proof. We expand $ and <3? _1 (|19p in the formal Laurent series 

$ L = a + ai5 _1 + a 2 S~ 2 + . . . , = /3 + AS" 1 + /? 2 5~ 2 + . . . . 

The coefficients £ Fq can be found recursively: 

l ( V«1.0 \ o _ c-i / Vll l,l 



1,1 / \ "»"1,0 



"i^^fe-fe^f^V ft=5-i(^-^5-i(5£M), (37) 



Then ()34|) follows from (|32p and formulas (|35p and (|36p follow from (|30p and (|3ip for n = 1 and 
n = 2 respectively. ■ 
There is no obstacle to compute equations corresponding to n = 3,4,5, etc. Of course, the 
expressions get bigger. Obviously, if we consider the recursion operators D\ in the T direction, we 
immediately get : 

Proposition 2. // a recursion operator 9\ is represented by a first order formal Laurent T series 
$Hl = T\T + fo + r_iT -1 H , then 

(i) (S - l)(ln h) = (T- 1)T- X (in , (38) 

(ii) (5-l)(f ) = (T-l)T- 1 (f 1 F), (39) 

(iii) (S-l)(r-iS-\h) + tf + hS(r-i)) = (T-l)(a 2 ), (40) 

w/iere 

& 2 = S-^h F) {S-Hro) + r - S~ 2 (n F) } - (1 + S" 1 ) (n (n f) 

cmc? F, G denote 

p_ Qux^T (Qwo.i ) ~~ Quq,qT (Qu lt i) g~ ^ u °'° 
Quq^T ^{Qui,i) Qu 01 
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It is not difficult to show that the existence of a formal series satisfying equation (|17p is the 
necessary condition for the existence of a hierarchy of higher symmetries (of increasing order) for a 
difference equation. Thus, formulas (|34H36p and (|38H40p can be seen as integrability conditions for 
the difference equation defined by function Q. Indeed, the right hand sides of these formulas are 
expressed in the terms of the difference equation only, and it is a rather non-trivial fact that the 
corresponding left hand sides are in the image space of the operators T — 1 or S — 1 respectively. 
If this is the case we can find the coefficients r\ , f i and proceed to the next conditions (|35p and 
(j39|) and so on. 

Theorem H] provides infinitely many integrability condition of that type. Integrability conditions 
(|34p and (|38|) have been earlier found by Levi and Yamilov [25] from the symmetry analysis of 
difference equations. They actually found two more conditions of this type, which in our setup 
correspond to the Taylor formal series (see footnote [1]) and can be easily recovered. Using the 
Taylor formal series we can double the number of the conditions, which would be appropriate for 
asymmetric or linearisable equations. 

It is not a straightforward exercise to apply the integrability conditions obtained. Nevertheless, 
they can be perfectly used for testing integrability as well as for classification of integrable difference 
equations (see for example [25]). 

4 Recursion operator, symmetries and conservation laws for the 
Viallet equation 

In this section we apply the theory developed in the previous sections to the case of the Viallet 
equation, as well as to the ABS equations. The first is characterized as integrable because its 
algebraic entropy vanishes, [20j, while the multidimensional consistency of latter implies their 
integrability, [17J. Here, we will show that they are integrable also in the sense that they admit 
infinite hierarchies of symmetries by constructing recursion operators for all of them. In order to 
make our presentation self-contained, we first introduce some properties of these equations. Next 
we present two symmetries and a recursion operator for the Viallet equation and demonstrate how 
it is related to corresponding operators for all of the ABS equations through an explicit example. 
Finally, we present some higher order conservation laws for all of the equations under consideration. 

Throughout this section we assume that, the defining function Q of equation ([1]) is irreducible. 
This means that parameters a, cannot take values for which Q can be factorized as a product of 
two polynomials. By making appropriate choices of the coefficients a% we can reduce the Viallet 
equation to every equation from the ABS list (see Appendix). The ABS equations depend on two 
so-called lattice parameters a and f3. 

Let us introduce two polynomials in terms of Q for equation (fTJ), which will be used in the rest 
of this section. They are determined by the discriminants of Q, 




(41) 
(42) 



Polynomials h(uofi, u^o), /i(uo.o> u o,i) are symmetric and biquadratic [18J. In the case of the ABS 
equations the polynomial h can be factorised and presented in the form 




(43) 
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where A; is a skew-symmetric function of the lattice parameters and / is symmetric and biquadratic 
polynomial of u and x [T7j. Polynomials f(u, x; a) for the ABS equations are listed in the Appendix. 

The covariance of these equations is another useful property for our presentation. For the 
Viallet equation this means that it is invariant under interchanging u±q and «o,i- F° r all of the 
ABS equations, covariance means that we have to interchange not only u\ t Q and Uo,i> but the lattice 
parameters a and /3 as well. This property requires us to study symmetries, recursion operators 
and conservation laws only for the one direction of the lattice (and recover the complementary set 
using the covariance). We shall assume that symmetries and canonical densities are elements of 
J- s and the recursion operator is a s-pseudo-difference operator. In most cases functions shifted 
by S only will contribute to the expressions. For this reason from now on, we will use one-index 
notation for the shifts of the polynomial h, i.e. 

h = h(u 0i0 ,Ui,o) , ^ = S l h(u 0fi ,ui fi ) , 

and symmetries, e.g. = S^(K^). We shall omit the index zero for unshifted functions. 

We now turn our attention to the symmetry analysis of the equations. It is known [21, 26j that 
these equations admit a generalised symmetry of order (—1, 1), which is given by the formula 

K W := h _ l du h = *z! + \d u _ w h^ . (44) 

Using Definition [5] and formula ([5]), we can directly verify the following result: 
Proposition 3. The Viallet equation (0) possesses a symmetry of order (—2, 2) which has the form 

= h J}zl ( I + I ) . (45) 

(U 1)0 - «-l,o) 2 \ "2,0 ~ "0,0 1*0,0 - "-2,0/ 

Actually, the Viallet equation possesses infinitely many local generalised symmetries. These 
hierarchies of symmetries can be constructed by applying successively a recursion operator on 
the seed or root symmetries and K^ 2 \ the starting points for a hierarchy of symmetries 



Theorem 5. The Viallet equation FT]) possesses a recursion operator 9\ = H o X, where operators 
H and X are defined as 

X - Is (46) 
h-xhhx s _ s _ x h_ihhi 



(«1,0 - M_l > o) 2 (u 2 ,0 - -U ,o) 2 ' ("1,0 - U-lfi) 2 {u2fl ~ UQfl) 2 

+ 2 S (S - l)- 1 o + 2 (S - l)- 1 o . (47) 

In this Theorem T~L and X are compatible Hamiltonian and symplectic operators for the Viallet 
equation, they will be discussed in details in our paper |30] including a complete proof of the 
Theorem. In the same paper, we will also show that both m e {K^) and ^ e (K^) are in T s for all 
<£N. This is not obvious at all, since the Hamiltonian operator 7i contains a pseudo-difference 
operator (S — 1) _1 . Symmetries and are called the seed (or root) symmetries for the 
Viallet equation. Starting from them, two infinite hierarchies of symmetries can be generated by 
applying the operator SH . Later in this Section we shall demonstrate the proof of Theorem [5] for 
the HI equation (Theorem [6]) . 

As an application, we now compute explicitly the action of $H on the first symmetry ifW . 
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Proposition 4. The Viallet equation (0) possesses a symmetry of order (—3,3) which is given by 



(ui j0 - n_i >0 ) 2 y («2,o - ^o,o) 2 (1*0,0 - w-2,0) 2 y 

+ ( + ^ . (48) 

V u 2,0 ~ uo,o u 0fi - n_ 2 ,o / 

Proof. To shorten the expression, we introduce the notation: 

w := , w k =S k (w). 

First we work out the product of operator % and X 

VK = hh- 1 w 2 {wfS 2 + S- 2 w 2 ) + 2K^KW His + S' 1 ^) 
-w 2 [h-! h lW \ + fc_ 2 hw 2 _ x ) + [k^K {2 \ + K®Kty 

+ 2^W(5-l)- 1 o (±=L +2^(5 -1)^0 £z± . (49) 

Acting with operator $H on the symmetry and using the identities 
K®ZK® = (5-1) (-^- K® K®) , i = 1, 2 ; 

K (2) XK (1) = (5 _ 1) fzl K (2)^(D + I(' 1 + 1 W) 

\« 2\ v n 2i o--u ,o Wo,o-M-2,o/ 
+ (5 - + 1) (- ^ ? 

\2(u 1)0 - u- 1)Q y(u ,o - u- 2 ,or 
K m XK w = -k<®ikV> + (5-1) f ±- kWk® + K^K (1 l); 

\h-i h-i ) 

we obtain expression (|48|) for symmetry i^ 3 ). ■ 
Using the above identities, one can easily compute the symmetry = D\(K^) of order 

(-4,4). Applying 9K N to the seeds and K^> we can generate two hierarchies of symmetries 
of order (-2JV - 1, 2N + 1) and (-27V - 2, 2iV + 2) respectively for the Viallet equation. 

Remark 1. Symmetries and recursion operators for the ABS equations can be obtained from the 
symmetries of the Viallet equation and the operator D\ = HI, given in Theorem by replacing 
the polynomial h ( everywhere in the symmetries and the operators I, H) with the corresponding 
polynomials f := f(uo 0, ui 0, a). For all ABS equations the polynomials f are listed in the Ap- 
pendix. Equations H1-H3, Ql and Q3g = o admit point symmetries ]2b^ . For these equations, we 
find additional Hamiltonian operators of zero Laurent order which are given in the Appendix. 

Proposition 5. If the recursion operator 9\, defined in Theorem^ is restricted to the HI equation 

(it ,o - 1*1,1) (1*1,0 - "0,1) = a - (3, (50) 
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then it becomes 
where 



Km = (5 + 1)(S - l)- 1 (5 - S- 1 ) . (51) 

"1,0 ""-1,0 "1,0 ~ "-1,0 

Proof. The restriction to the HI equation implies that h = 1. Thus, the first two symmetries 
and (|45p become 

= and = -, ^ ( + 



"1,0 -"-1,0 ("1,0 - u -l,o) 2 \ u 2,0 ~ "0,0 "0,0 -"-2,0/ 

respectively. Substituting them into (|49p . a direct calculation leads to the statement. ■ 
The operator given by ()5ip is a recursion operator of the HI equation. To prove it, we will use 
the following result. 

Lemma 1. Operator 

11 ° 

= -, + 



("2,1 -"0,l) 2 ("1,0 - "-l,o) 2 ("1,1 - "-l,l)("2,l - "0,1 ) 

+(q - /3) ' ' : 



- u_i i i) 2 (ni j0 - n ,i) 2 (it 2 ,i - "o,i) 2 ("2,o - "1,1) 

1 a-fl 
+ 



"1,0 - "-1,0 ("2,0 - "0,o)("l,0 - "0,l) 

(S-l)-lo(^ , (52) 

V"2,0 - "0,0 ("1,0 - "-l,0j("0,0 - "-1,1 J / 

satisfies the following identity for all solutions of the HI equation {5U\): 



1 + 7 \ 2 S ) °*Hi = JWofl + - - \ 2 S); (53) 
(iti,o-"o,i) / V ("i,o-"o,i) 



S + ( \ 2 ) ° T(m m ) = Wlo(s+ a P ) . (54) 
("1,0 - "o,i)V V ("1,0 - "o,i)v 

Proof. Writing out explicitly both sides of equation (|53p . we factorise terms in both sides 
according to the different powers of S. Collecting the coefficients of the different powers of S in 
both sides of the resulting expression, their equality follows by putting them into J-q. In the same 
way, we can prove identity ([Ml)- B 



Theorem 6. Operator VKhi given by \51\) is a recursion operator of the HI equation S50\) with 

V = ("1,0 - "0,1) ° ("i,o - "o,i) _1 • (55) 

Proof. From Theorem [1] it follows that, 9\hi is a recursion operator of the HI equation 
provided both formulas (|17p and (|18p are valid for operator We rewrite them as 

(Q uifi s + Q« ,o) ° 9tin = (Qu lfi s + Q« ,o); (56) 

(<9 U i.i«S + Q U0 . 1 )oT(5Hhi) =¥°(Qu 1 , 1 S + Q U01 ). (57) 
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For the HI equation, we have 



a. — P 

Qui o S + Qu , = («o,o - ui,i) S + (ni j0 - n ,i) = (ni j0 - uq,i) ( 1 + 



Qu x i 5 + Q no j = - (ui,o - ^0,1) S - (n ,o - = - («i,o - ""0,1) 5 + 



(wi j0 - n ,i) 2 
a — /3 



("1,0 - "0,1 



>2 



Using Lemma [U we can easily obtain (|56p and ()57p and thus we proved the statement. ■ 
Having constructed a recursion operator for equation ([1]), one can use the results of Section 
13.21 to find conserved densities for this equation. If we apply directly Theorem 2] to the recursion 
operator (|49p . since the Laurent order of the recursion operator 91 is two, we will get only the 
"even" part of the infinite sequence of the canonical conservation laws . Although it is impossible 
to compute a square root of a generic difference formal series, for the recursion operator (|49p it can 
be done and thus one can find a first order formal recursion operator 

*= , h ~ l , 2 S+ 2Kil) -d U0)0 KW + ... , 

(Ul,0-U-1,0) "0,0 -"-2,0 

such that 9\ 2 = 9\. The first order recursion operator provides us with a complete set of canonical 
conservation laws and the first three densities of this sequence are 

Po = In (<9 Uli0 K (1) ) = hi/i-i - 2 ln(n lj0 - u_i, ) , (58) 



Pi 



"0,0 ~ "-2,0 



<o^ (1) , (59) 



P, = . ^ ,A ( k - lhl v + ( k - 2k , , )+^(KM K V + KWK%. m 
[uifl - u-i,or \{u 2 ,o - "o,or ("0,0 - u-2,0) J h-i \ ) 

The first and the second densities are the logarithmic residue and the residue of 91, respectively, 
while P2 is the residue of the recursion operator = S^l 2 (|49p . 
According to Proposition [H the corresponding fluxes are 



00 



0" 2 



W^tK), (61) 

^(r^^ (62) 
^ (m,o - «-i,o) / 

0-1 / 2F s _ 2 ( Fhh-x 



where 



h J \("2,o - "o,o) 2 ("i,o - "-1,0) 
-(1 + O ^ , (63) 

VQiti.o ("1,0 - "-i,or("o,o - u-2,0) J 



and Q is the defining polynomial of equation ([!]) (or the ABS equations). 
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The covariance of the equation allows us to produce three more conservation laws simply by 
changing u^o to Uq^, h to h and operator S to T in the above densities and fluxes. On the other 
hand, these conservation laws can also be restricted to the ABS equations in the way described in 
Remark [TJ 

For the HI equation, its recursion operator (|5ip is first order and we can use relations (|33p to 
derive some conserved densities. Actually, the following formulas for conserved densities are valid 
for the equations admitting Hamiltonian operators of zero Laurent order, i.e. equations H1-H3, QI 
and Q3,5 = o. The conserved densities for all these equations can be written in a relatevly compact 
form (they are equivalent to the ones given in Proposition [1]): 

ln(5- 1 /)-21n(n 1>0 -n_ 1 .o), 

/ 

(Ulfi - u- lfi )(u 2 ,o - 160,0)' 
(0i + S( 0l )) 2 -S( ei ) 2 , 
(Qi + S(ei)) 3 - S( ei ) 3 + 3^5(^i)5 2 ( ft ), 

(01 + <s(ei)) 4 - s( ei ) 4 + 4 e i5( e i)5 2 ( ei )(ei + 25(ei) + s\ 01 ) + s 3 ( 01 )). 

For the HI equation (|50p . some of the above densities were found in [23J by the Gardner method. 

5 Conclusions 

The recursion operators presented in this paper produce infinite hierarchies of symmetries and 
conservation laws for all ABS equations. For the HI equation we have shown that the proposed 
recursion operator, 9\hi , does satisfy defining equation (fT6j) and thus acting on its symmetry 
it will produce its new symmetry 9^h 1 [K^ 1 '), provided such action is well defined (i.e. the result 
belongs to Tq). The proof that we can repeatedly apply the recursion operator and each time get 
a symmetry requires us to investigate the property of Nijenhuis operators [311 132] and to develop 
a construction similar to the famous Lenard scheme for the Korteweg-de Vries equation [33, 34J. 
The generalisation of Lenard's scheme for difference equation as well as a complete proof that the 
operator 9\ defined in the Theorem [5] does satisfy (|16p for the Viallet equation will be published in 
our paper [30] . 

For a quadrilateral equation the first integrability condition has been found by Levi and Yamilov 
[25] • In this paper we have proposed an infinite sequence of integrability conditions in the form of 
canonical conservation laws. The Levi- Yamilov condition is the first element of this sequence. We 
have shown that canonical conservation laws are necessary conditions for the existence of a formal 
recursion operator for a difference equation. As in the case of partial differential equations, cf. 
[3j El [6], one can prove that the existence of a formal recursion operator follows from the existence 
of an infinite hierarchy of symmetries of increasing order for a difference equation (a general and 
detailed proof of this assertion will be published elsewhere) . There are no doubts that integrability 
conditions can be successfully used for testing and classification of integrable difference equations 
(see for example |25J). 

Our theory can be extended to difference equations of "higher order", i.e. equations on a 
rectangular with function Q depending on the set of variables {u Ptq |0 < p < N , < q < M}. 
For the existence of the elimination map £ we should assume that equation Q = can be uniquely 
resolved with respect to the corner variables ito,0) ^iV.O) uo,M and un,m- In this case the dynamical 



So = 

Qi = 

Q2 = 

Q3 = 

Qa = 
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variables can be denned as 



U s = {u n;P \neZ, 0<p<M}, U t = {u qiTn \meZ, 0<q<N}, U = U s UU t . 

We are working on the extension of our theory to the higer order case, the vector case, as well as 
to the case when Q is not assumed to be an afhne-linear polynomial. 
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Appendix : The ABS equations 

For every equation in the ABS classification [T7] we present its defining polynomial Q, the corre- 
sponding choices for in ([1]), the corresponding polynomial / and factor k (|43p . the Hamiltonian 
operator H and the seed (root) symmetry of a hierarchy. 

The recursion operator for all of the ABS equations can be written as 9\ = % o X, where the 
symplectic operator Z is 

and the Hamiltonian operator H is either 

^ f-i f fi ^ _ £-1 /-l / fi 



(Ul,0 ~ U-lfi) 2 (u2,0 ~ UOfl) 2 ' (Ulfl ~ U-i )0 ) 2 (u 2 fl - U fl) 2 

+ 2 ( KM S(S- l)- 1 ^ + {S - l)~ l K^ ) , (64) 



which is valid for all of the ABS equations, or the one given in the list for the case of H1-H3, 
Ql and Q3,5 = o. In (|64p and in the lists, and denote the first symmetries of the ABS 

equations, which have the following form 

K {1) = f - \d Ul J = ^ + 

"1,0 - M-1,0 2 "1,0 - U-ifi 2 



("1,0 - «-l,o) \U2,0-U ,0 U 0i0 -U-2,0 



HI 

Equation Q = (uo,o — ui,i)(u\fi — uo,i) — a + f3 

Parameters a% = —a^ = 1, 07 = — a, a\ = 02 = 04 = «6 = 

Polynomial h f = 1 

k = f3 — a 

U U = KW(S + 1)(S - l)~ l KV-) 
Seed symmetry 
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H2 

Equation Q = (u ,o - i*l,i)(«l,o - ^0,1) + W ~ «)(«o,o + «l,o + ""0,1 + - a 2 + /3 2 

Parameters 03 = —05 = 1, a§ = j3 — a, 07 = /3 2 — a 2 , ai = 02 = CI4 = 

Polynomial ft / = 2(ito,o + ^1,0 + <*) 
k = f3 — a 

n H = K<V(S + 1)(S - i)- 1 ^ 1 ) - (S - 1)(S + 

Seed symmetry 
H3 

Equation Q = a(uo,oUi,o + ^0,1^1,1) ~~ /S(«o,o«o,l + "1, o«l,l) + ^(° 2 ~~ /^ 2 ) 

Parameters 03 = a, 05 = — /3, 07 = 5(a 2 — /3 2 ), a\ = CL2 = = = 

Polynomial ft f = uq oui + 016 
k = a 2 '-p 2 

Seed symmetry 

Qi 

Equation Q = a(u ,o - ^o,i)( n i,o - u i,i) ~ /3(uo,o ~ «i,o)(«o,l - ui,i) + S 2 a/3(a - 0) 

Parameters (13 = a, 04 = /3 — a, 05 = — /3, 07 = <5 2 a/3(a — /3) 

Polynomial ft f = ((«o,o — ""l^) 2 — a 2 5 2 )/a 
k = a(3(/3 — a) 

H n = KW(S + 1)(S - l)" 1 ^ 1 ) - 5 2 (S + 1)(S - l)- 1 
Seed symmetry 

Q2 

Equation Q = a(u 0fi - «o,i)( n i,o - «i,i) - /3(«o,o - ^i,o)Oo,i - «i,i) 

+a/3(a - /0)(uo,o + ^1,0 + ""0,1 + «i,i) ~~ - /3)(a 2 - a/3 + /3 2 ) 

Parameters 03 = a, 04 = /3 — a, 05 = — /3, a6 = a/3(a — /3), 07 = a/3(/3 — a)(a 2 — a/3 + /3 2 ) 

Polynomial ft / = ((u ,o - ui,o) 2 - 2a 2 (u ,o + 1*1,0) + a 4 )/a 
/c = a/3(/3 — a) 
given by (fM|) 

Seed symmetries i^ 1 ), if( 2 ) 

Q3 

Equation Q = (/3 2 - a 2 )(« ,o«l,l + Ml,o«o,l) + ~ l)(«o,o«l,o + «0,iwi,l) 

iu , \ 5 2 (a 2 -/3 2 )(a 2 -l)(/3 2 -l) 
-a(/3 z - l)(lt ,0«0,l + «l,0«l,l) " 4a/3 i 

Parameters a 3 = /3(a 2 - 1), a 4 = /3 2 - a 2 , a 5 = -a(/3 2 - 1), 07 = _ W^gzlK^li 

Polynomial ft / = 4a ^i_^ (4a(«n ,o - iti i0 )(mti,o - uq,o) - (a 2 - 1) 2 5 2 ) 

k = (a 2 -/3 2 )(a 2 -l)(/3 2 -l) 
H (i.) 5 = 0:^ = + - l)" 1 ^ 1 ) - in 0i0 (5 + 1)(«S - l)~Vo 

(ii.) 5/0: given by 
Seed symmetries (i.) 5 = 0: ifW; (ii.) 5/0: K^, 
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Q4 

Equation Q = sn(a)(uo,oni,o + ito,i«i,l) ~ sn(/3)(u 0) o«o,l + wi,o«i,l) 

-sn(a - /0)(«o,owi,i + «l,o"o,i) +ksn(a) sn(/3) sn(a - 0) (1 + ^0,0^1,0^0,1^1,1) 

Parameters a\ = aj = k sn(a) sn(/3) sn(a — 03 = sn(a), 04 = — sn(a — 05 = — sn(/3) 

Polynomial h f = -ksn(a)(l + ng nf ) + i5 ^y (u§ )0 + u\ fl - 2cn(a)dn(a)u ,o«i,o) 
k = sn(a)sn(/3)sn(/3 — a) 

7^ given by (O 

Seed symmetries X^ 1 -', 
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